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Abstract. In this paper we generalize the involutive methods and algorithms de- 
vised for polynomial ideals to differential ones generated by a finite set of linear 
differential polynomials in the differential polynomial ring over a zero characteristic 
differential field. Given a ranking of derivative terms and an involutive division, we 
formulate the involutivity conditions which form a basis of involutive algorithms. 
We present an algorithm for computation of a minimal involutive differential basis. 
Its correctness and termination hold for any constructive and noetherian involutive 
division. As two important applications we consider posing of an initial value prob- 
lem for a linear differential system providing uniqueness of its solution and the Lie 
symmetry analysis of nonlinear differential equations. In particular, this allows to 
determine the structure of arbitrariness in general solution of linear systems and 
thereby to find the size of symmetry group. 



1 Introduction 

Among the properties of systems of analytical partial differential equations 
(PDEs) which may be investigated without their explicit integration there are 
compatibility and formulation of an initial-value problem providing existence 
and uniqueness of the solution. The classical Cauchy-Kowalevsky theorem 
establishes a certain class of quasilinear PDEs which admit posing such an 
initial- value problem. The main obstacle in investigating other classes of PDE 
systems of some given order q is existence of integrability conditions, that is, 
such relations for derivatives of order < q which are differential but not pure 
algebraic consequences of equations in the system. 

An involutive system of PDEs has all the integrability conditions incor- 
porated in it. This means that prolongations of the system do not reveal 
integrability conditions. Extension of a system by its integrability conditions 
is called completion. The concept of involutivity was invented hundred years 
ago by E.Cartan |l| in his investigation of the Pfaff type equations in to- 
tal differentials. For these purposes he used the exterior calculus developed 
by himself. The Cartan approach was generalized by KahlerQ to arbitrary 
systems of exterior differential equations. The underlying completion proce- 
dure 1^ was implemented in 
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In his study of the formal power series solutions of PDEs, Riquier intro- 
duced Q a class of relevant rankings for partial derivatives and considered 
systems of orthonomic equations which are solved with respect to the highest 
rank derivatives called principal. Thereby, these derivatives, by the equations 
in the system, are defined in terms of the other derivatives called parametric. 
An integrability condition gives a constraint for parametric derivatives, and 
that of them of the highest ranking becomes the principal derivative. Recently 
Riquier 's class of rankings was generalized in |Q. 

Janet made the further development of Riquier's approach. He observed 

that the integrability conditions may occur only from prolongations with 
respect to certain independent variables called nonmultiplicative. Prolonga- 
tions with respect to the rest of variables called multiplicative never lead to 
integrability conditions. Given a set of principal derivatives, Janet gave the 
prescription how to separate variables into multiplicative and nonmultiplica- 
tive for every equation in the system. He formulated, on this ground, the 
involutivity conditions for orthonomic systems and designed an algorithm for 
their completion. This approach to completion is known as Riquier-Janet 
theory and was implemented in [p| pO| , pT[ |. 

A system satisfying the Janet involutivity conditions is often called pas- 
sive. This involutivity is generally coordinate dependent. On the other hand, 
the modern formal theory of PDEs developed in 60s-70s by Spencer and 
others (see [ p^HT^ ) allows to formulate the involutivity intrinsically, in a co- 
ordinate independent way. The formal theory relies on another definition of 
multiplicative and nonmultiplicative variables which was known to Janet as 
long ago as in 20s, but called nowadays after Pommaret because of its im- 
portance in the technique presented in . The implementation in Axiom of 
completion based on the formal theory was presented in . 

Thomas in [^j used another separation of independent variables into mul- 
tiplicative and nonmultiplicative and generalized the Riquier-Janet theory to 
non-orthonomic algebraic PDEs. Given a system of PDEs, he showed that in 
a finite number of steps one can: (i) check its compatibility; (ii) if the system 
is compatible, then split it into a finite number of simple systems involving 
generally both equations and inequalities and such that their equation parts 
are orthonomic and can be completed to involution. This splitting is similar 
to that generated by the Rosenfeld-Grobner algorithm JlTf . 

In paper ||l8|] for Pommaret separation of independent variables it was 
shown that involutive (passive) basis of a non-differential polynomial ideal is 
a Grobner basis. The implementation in Reduce of the proposed completion 
algorithm for polynomial bases demonstrated a high computational efficiency 
of the involutive technique. However, Pommaret bases may not exist for pos- 
itive dimensional ideals unlike Janet and Thomas bases. 

The above classical separations of variables into multiplicative and non- 
multiplicative are particular cases of involutive monomial division, a concept 
invented and analyzed in jl^ (cf. ||2^). The polynomial completion algo- 
rithms designed for a general involutive division ||l^ , pT| were implemented in 
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Reduce for Pommaret division. Different involutive divisions and completion 
of monomial sets have also been implemented in Mathematica ||2^. In 
we generalized the algorithm of paper pl| to arbitrary completion ordering. 

One more efficient method for the completion of linear PDEs to an involu- 
tive form called standard which is not based on the separation of variables was 
developed in |Q and implemented in Maple. The extension of this method 
to nonlinear PDEs is given in paper |25[ . 

In the present paper we generalize the involutive methods and algorithms 
devised in 19, 2l] , p3| for polynomial ideals to differential ideals generated by 



a finite set of linear polynomials. We formulate the involutivity conditions for 
the differential case. If a set satisfies the involutivity conditions it is called 
an involutive basis. Similar to the pure algebraic case, a linear involutive 
basis is a differential Grobner basis [^6[|27|] which is not generally reduced. 
We present an algorithm for computation of a minimal involutive basis. This 
algorithm is the straightforward generalization of the polynomial involutive 
algorithm |2^,^. As well as for the latter, the correctness and termination 
of the former hold for any constructive and noetherian involutive division. 

An important application of the involutive method is posing an initial 
value problem providing the unique solution of a system of PDEs. For linear 
involutive systems we formulate such an initial value problem and thereby 
generalize the classical results of Janet [^j to arbitrary involutive divisions. 
This formulation makes it possible, among other things, to reveal the struc- 
ture of arbitrariness in general solution. Given a linear involutive basis, we 
write also the explicit formulae for the Hilbert function and the Hilbert poly- 
nomial of the corresponding differential ideal which are the straightforward 
generalizations of their polynomial analogues ||2^,^ . 

Another important application of the new algorithm is the Lie symmetry 
analysis of nonlinear differential equations. It is because of the fact that com- 
pletion to involution is the most general and universal method of integrating 
the determining system of linear PDEs for infinitesimal Lie symmetry gen- 
erators [^ . Moreover, an involutive form of determining equations allows to 
construct the Lie symmetry algebra without their explicit integration |^ . In 
particular, for an involutive determining system the size of symmetry group 
can easily be found that was shown for Janet bases in . Though reduced 
Grobner bases for the determining equations do not generally reveal infor- 
mation on Lie symmetry groups, and more generally on the solution space, 
so explicitly as involutive bases, they are also very useful for Lie symmetry 
analysis as shown in . The facilities of the Maple package devised by the 
first author and used in the paper go far beyond linear differential systems, 
and it can also be fruitfully applied to nonlinear systems. 

2 Preliminaries 

Let M = . . . , ym} be a differential polynomial ring [|3iy32|] with the set 

of differential indeterminates {yi, . . . , j/m}, and IK C M is a differential field 
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of zero characteristic with a finite number of mutually commuting derivation 
operators . . . , 9/9a;„. Elements in M are differential polynomials in 

{j/i: ■ ■ • I Vm}- In this paper we use the following notations and conventions: 

f,g,h,p(z M are linear differential polynomials. 

F,G, H cM. are finite sets of linear differential polynomials. 
= {/ = I / e i^} is a linear system of PDEs. 

N={0,f,2,...}is the set of nonnegative integers. 

a,/3, 7 £ N" are multiindices. 

Icm^a, [3) is the least common multiple of alpha, beta. 

X = {xi, . . . , Xn} is the set of independent variables. 

R ~ K[X] is the polynomial ring over the field K of zero characteristic. 

R^M = {x" = xl"- ■ ■ ■x'^^ I ai e N} is the set of monomials in X. 

i — f , . . . , n indexes derivation operators di = d/dxi- 

j — 1, . . . ,m indexes indeterminates yj . 

u, V, w are elements in M. 

?7, C M are finite monomial sets. 

[U] is the monomial ideal in R generated by U . 

degi{u) is the degree of Xi in w G M. 

deg(u) = X)"=i degi(u) is the total degree of u. 

^^yj = dli^-ali^ yj ^ derivative. 

ord^daUj) = '^i is the order of dayj- 

9, ■& are derivatives. 

Icm{dayj,di3yj) = dicm(a,l3)yj- 

^, -<c are rankings of derivatives. 

ld{f) is the leading derivative in / e M. 

lc{f) G if is the coefficient of ld{f). 

ld{F) is the set of leading derivatives in C M. 

[F] is the differential ideal in M generated by F. 

L is an involutive division. 

L{u, U) is the set of (L— )multiplicative monomials for u E U. 

NMl{u, U) is the set of L— nonmultiplicative variables for u £ U. 

Ml{u, U) is the set of L— multiplicative variables for u £ U. 

a;" £ M is the monomial associated with the derivative dayj. 

UjLiUj is the monomial set associated with the set ld{F) = U"Li{ldj{F)}. 

W = U^Li{W^j I Wj C M} is the complementary set for U'J'^iUj. 

G is the set of L— generators of W. 

•d — OlO is a multiplicative prolongation of 9. 

■ 9 is the nonmultiplicative prolongation of 9 w.r.t. Xi. 
da • 9 is a nonmultiplicative prolongation of 9. 
NFl{p, F) is the L— normal form of p modulo F. 

NML{f, F) C X is the set of (L— )nonmultiplicative variables for f E F. 

Cl{F) — 'J9£id{F){^ I = dL9} is the cone generated by F. 
In this paper we distinguish two rankings (c.f. p3| ): a main ranking and a 
completion ranking denoted by >~ and >-c, respectively. The main ranking 
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will be used, as usually, for isolation of the leading derivatives in differen- 
tial polynomials whereas the completion ranking serves for taking the lowest 
nonmultiplicative prolongations by the normal strategy [l9[ | and thereby con- 
trolling the property of partial involutivity introduced in Sect. 4. 

3 Basic Concepts and Definitions 

Throughout this paper we exploit the well-known algorithmic similarities 
between pure algebraic polynomial systems and linear differential systems 
[ P^IJb^ . In so doing, the basic algorithmic ideas go back to Janet j|] who 
invented the constructive approach to study of PDEs in terms of the corre- 
sponding monomial sets which is based on the following association between 
derivatives and monomials: 

aaiH hQ,i,, . 

~ dx""^ ■ ■ ■ dx"'" X — x^ x^ . yj^) 

The monomials associated with the different indeterminates yj are to be con- 
sidered as belonging to different monomial sets Uj G M indexed by subscript 
j of the indeterminate. 

Definition 1. ||3^ A total ordering -< over the set of derivatives daUj is 
called a ranking if it satisfies: (i) didaUj >- daUj, (ii) daTJj >~ dpyk 
d-fdaVj >- dydpyk for all i, j, k, a, (3,"f. A ranking -< is said to be orderly if 
9 >- d whenever ord{9) > ord^d). 

The association (|^) implies the reduction of a ranking -< to the associated 
admissible monomial ordering, and throughout the paper we shall assume 
that 

di ^ ■ ■ ■ >~ dn <J=> Xi)^ X2> ■ ■ ■ >~ Xn ■ (2) 

Remark 2. Given a finite set F C M and a ranking set ld{F) of the leading 
derivatives is partitioned ld{F) = [Jjldj{F) into subsets ldj{F) corresponding 
to different indeterminates yj which occur in ld{F). For an involutive division 
L defined as follows each subset generates for every its element the separation 
of independent variables into multiplicative and nonmultiplicative ones. 



Definition 3. [19 An involutive division L on M is given, if for any finite 
monomial set f/ C M and for any u € U there is given a submonoid L(u, U) 
of M satisfying the conditions: 

(a) If G L{u, U) and v\w, then v G L{u, U). 

(b) If u, w G [/ and uL{u, U) n vL{v, U) ^ 0, then u G vL{v, U) 
or w G uL{u, U). 

(c) If u G t/ and w G uL{u, U), then L{v, U) C L{u, U). 

(d) If F C t/, then L(ii, U) C L{u, V) for aU w G F. 
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Elements of L{u, U) are called multiplicative for u.liw ^ uL{u, [/), m is called 
an involutive divisor or (L—) divisor of w. In such an event the monomial 
V = w/u is called L— multiplicative for u. If m is a conventional divisor of w 
but not L— divisor, then v is called nonmultiplicative for u. 

Remark 4- Definition ^ for every m G [/ provides the partition 

X ^ Ml{u,U)UNMl{u,U), MlCL{u,U) (3) 

of the set of variables X — {xi, . . . , Xn} into subset Ml{u, U) of multiplica- 
tive variables for u and subset NMj^iu, U) of the remaining nonmultiplicative 
variables. Conversely, if for any finite set J7 C M and any u U the parti- 
tion of variables into multiplicative and nonmultiplicative is given such 
that the corresponding submonoid L{u) satisfies the conditions (b)-(d) in 
Definition ^, then the partition generates an involutive division. 

Definition 5. A monomial set U is called L—autoreduced if uL{u, U) fl 
vL{v, [/) = holds for all distinct u,v £U . 

Definition 6. ||l^ A monomial set U is called an L— completion of a set 

u (zjj a 

(Vii e U) (Vw e M) {3v eU) [uwe vL{v, U) ] . 

If there exists a finite L— completion [/ of a finite set U, then the latter 
is called finitely generated with respect to L. The involutive division L is 
noetherian if every finite set U is finitely generated with respect to L. If 
U — U, then U is called L— complete. An L—autoreduced and complete set 
is called {L—)involutive. 

Definition 7. Given a monomial set U, the set Uueu m M is called the 
cone generated by U and denoted by C{U). The set U„g[/ U) is called 

the involutive cone of U with respect to L and denoted by Cl{U). 

Thus, the set U is an L-completion of U if C(J7) = CL(i/) = C(C/). Corre- 
spondingly, for an involutive set U the equality C{U) = Cl{U) holds. 

Whereas noetherity provides existence of a finite involutive basis for any 
polynomial ideal, another important properties of an involutive division called 
continuity and constructivity provide the algorithmic construction of involu- 
tive bases [ p^ . Continuity implies involutivity when the local involutivity 
holds whereas constructivity strengthens continuity and allows to compute 
involutive bases by sequential examination of single nonmultiplicative pro- 
longations only. We refer to papers |ll9| , pTlj2^ for description of these topics 
in detail. In those papers some examples of involutive divisions were studied 
(see also [Q) which include three divisions called after Janet, Thomas and 
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Pommaret, because they have used the corresponding separations of vari- 



p| Jl^ , p^ . Other two divisions called 
and a class of involutive divisions 



ables for involutivity analysis of PDEs 
Division I and II were introduced in ||^ 
called Induced division, since every division in the class is induced by an ad- 
missible monomial orderings, was introduced in p3[ . All those divisions are 
constructive and, except Pommaret division, they are noetherian. Below we 
use three of those divisions defined as follows. 

Definition 8. Janet division Let U C M be a finite set. Divide U into 
groups labeled by non-negative integers ai, . . . , (1 < i < n) : 

[ai, . . . ,ai\ — { u ^ U \ Uj — degj{u), I < j < i }■ 

Then Xi is multiplicative for u G U if i ~ I and degi{u) ~ max{(iegi(u) | v G 
C/}, or It e [ai, . . . , ai-i] and degi{u) = max{(ie(7i(v) | v £ [ai, . . . , cti-i]} for 
i > 1. 



Definition 9. Pommaret division |]12|. For a monomial u = x"^ ■ ■ ■ x'^'' with 
afc > the variables Xj , j > k are considered as multiplicative and the other 
variables as nonmultiplicative. For w = 1 all the variables are multiplicative. 



Definition 10. Lexicographically induced division [g3|. A variable Xi is non- 
multiplicative for u e [/ if there is v E U such that v -<Lex u and degi{u) < 
degi{v), where )^Lex denotes the lexicographical ordering. 

In the sequel Janet, Pommaret and Lexicographically induced divisions will 
be distinguished by the subscripts J,P and Dl^x, respectively. 

Example 11. Separation of variables for set U = {a;^a;3, a;ia;2, a;ia;|} and or- 
dering (j^) for the above defined three divisions: 



Element 
in U 


Separation of variables 


Janet 


Pommaret 


Lex. induced 


Mj 


NMj 


Mp 


NMp 






x'ix'i 


Xi,X2,X3 




X3 


Xl,X2 


Xl 


X2,X3 


XlX2 


X2,X3 


Xl 


X2,X3 


Xl 


Xl,X2 


X3 


xixl 


X3 


Xi,X2 


X3 


Xi,X2 


Xi,X2,X3 





The corresponding L-completions of U are 

ilj = {x{x3,xiX2,xixl,xlx2}, 

TT < 2 2 2 i+2 j+2 1 

Up = {XiX3, XiX2, XiX^, XiX2, . .. ,X{^ X2,-. ■ , x{ X3, . . .}, 

UoLe^ = {x\x3,XiX2,Xixl,XiX2X3). 

where i,j G N. This example explicitly shows the non-noctherity of Pom- 
maret division. 



8 



Definition 12. Given a finite set F C M, a ranking >- and an involutive di- 
vision L, the derivative = dpyj will be called a multiplicative prolongation 
oi = daUj G ldj{F) and denoted by = OlO, if the associated monomials 
satisfy € Uj). Otherwise the prolongation will be called nonmul- 

tiplicative. Respectively, the corresponding prolongation dpf of the element 
f & F with ld{f) — daHj will be called multiplicative and denoted by dL^f) 
or nonmultiplicative. The set Cl{F) — ^eeid(F){'& I — QlO} will be called 
the L—cone generated by F. If dif is a nonmultiplicative prolongation of 
f e F,we shall write G NML{f, F). 

4 Linear Involutive Differential Bases 

In this section we generalize the results obtained in papers [|9|j2^ for com- 
mutative algebra to differential algebra of linear polynomials. Proofs of the 
theorems are omitted because of similarity with the proofs of their algebraic 
analogues. 

Definition 13. Given an involutive division L, a finite set F C M of linear 
differential polynomials, a ranking >- and a linear polynomial p e M, we shall 
say: 

1. p is L— reducible modulo f £ F if p has a term a9, [a E K \ {0}) such 
that = dLld{f ). It yields the L— reduction p ^ g = p — {a/lc{f ))dpf 
where dpld{f) = 9. 

2. p is L~reducible modulo F if there is f G F such that p is L— reducible 
modulo /. 

3. p is in L— normal form modulo F, if p is not L— reducible modulo F. 

We denote the L— normal form of p modulo F by NFl{p, F). 

As a i— normal form algorithm one can use the following differential ana- 
logue of the polynomial normal form algorithm ||l9|] : 

Algorithm InvolutiveNormalForm: 

Input: p, F, L, ^ 
Output: h = NFl{p,F) 
begin 

h := p 

while exist f E F and a term a 9 (a G IfC \ {0}) of h 

such that 9 = dLld{f) do 

ciioose the first such / 

h:^h- {a/lc{f)) dpf where d0ld{f) = 9 
end 
end 
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Correctness and termination of this algorithm is an obvious consequence of 
Definition 13 and correctness and termination of the polynomial L— normal 
form algorithm [ p^ . 

Definition 14. A finite set F is called L—autoreduced if every f G F is 
irreducible modulo any other element g e F. An L—autoreduced set F is 
called {L—)involutive if 

(V/ e F) (Va e N") [ NF^d^f, F) = ]. 

Given a derivative i? and an i— autoreduced set F, if there exist f ^ F such 
that ld{f ) -<c I? and 

(V/ e F) (Va e N") ^) [ NF^d^f, F) = 0], (4) 

then F is called partially involutive up to the derivative {} with respect to the 
ranking -<c. F is still said to be partially involutive up to d \i § -<c ld{f) for 
all f £F. 



Corollary 15. If F C ^ is an L— involutive set, then every monomial set 
Uj G M (1 < J < to) associated with ldj{F) is L~involutive. 

Proof. It follows immediately from Definitions ^ and |l^. □ 

Theorem 16. An L—autoreduced set _F C M is involutive with respect to a 
continuous involutive division L iff the following (local) involutivity condi- 
tions hold 

(V/ e F) (Vx, e NML{f. F)) [ NFL{d^^ • /, F) = ] . 
Correspondingly, partial involutivity holds iff 
(V/ e F) {\/x, e NMLif, F)) (9,, • IdU) <c i9) [ NFL{d^^ • /, F) = ] . 



Theorem 17. If F d is an L— involutive basis of [F], then it is also a 
differential Grobner basis. 

The following theorem and corollary give an involutive analogue of Buch- 



berger chain criterion |35 in application to linear differential bases. 



Theorem 18. Let F be a finite L—autoreduced set of linear differential poly- 
nomials with respect to a continuous involutive division L, and NF^lp, F) be 
an algorithm of L— normal form. Then the following are equivalent: 
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1. F is an L—involutive differential basis of [F]. 

2. For all g e F,x e NML{g,F) there is f E F satisfying ■ ld{g) = 
dhldif) cind a ehain of elements in F of the form 

f = fk, fk-l, ■ ■ ■ , fo, 90, ■ ■ ■ , gm-l,gm = 9 

such that 

NFl {SL{f^-ufi).F) = NFl iSifo,go),F) = NFl {SL{g,-ug,),F) = 

where <i < k and < j < m, S'(/o,5o) the conventional differential 
S-polynomial Jg^ and SL^fi, fj) = dx ■ fi — dhfj is its special form which 
occurs in involutive algorithms. 



Corollary 19. Let F be a finite L— autoreduced set, and let dx • g be a non- 
multiplicative prolongation of g G F . If the following holds 



[Mh E F) dc,) idjd{h)-u^cldig-x)) [ NFLih ■ u, F) ^ ] , 
(3/,/o,5oeF) 



Idif) = dfjldifo) , ld{g) = d^ld{go) 
dx ■ ld{g) = dLld{f) , Icm {ld{fo), ld{go)) <c dx ■ ld{g) 
_NFL{dp-fo,F)^NFL{d^-go,F)=0 

then the prolongation dx ■ g may be discarded in the course of an involutive 
algorithm. 



5 Completion Algorithm 



The below given algorithm MinimalLinearlnvolutiveBasis is a differen- 
tial analogue of the polynomial algorithm MinimallnvolutiveBasis of pa- 
per |2^. In so doing, the conventional (non- involutive) autoreduction which 
is performed in line 2 of the latter algorithm omitted, as this autoreduction 
is optional 

Validity of the involutive chain criterion used in lines 11 and 23 is provided 
by Theorem |l^ and Corollary |l^. The proof of correctness and termination 
of the differential algorithm is identical to the proof for its polynomial ana- 
logue pTp^ . It follows, that if the main ranking >~ is orderly, then, given a 
generating set of linear differential polynomials and a constructive involutive 
division, algorithm MinimalLinearlnvolutiveBasis computes a minimal 
differential basis whenever the latter exists. If the division is noetherian, the 
basis is computed for any main ranking. 

Though the output basis for a noetherian division does not depend on the 
completion ranking, the proper choice of the latter may increase efficiency of 
computation. 
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Remark 20. If the algorithm MinimalLinearlnvolutiveBasis takes a con- 
ventional differential Grobner basis of the ideal [F] as an input, then it pro- 
duces the minimal involutive differential basis just by enlargement of the 
input set with its irreducible nonmultiplicative prolongations if any. This 
enlargement is done in the lower while-loop. 

Algorithm MinimalLinearlnvolutiveBasis 



Input: F, L, >- (main ranking), >-c (completion ranking) 
Output: G, a minimal involutive basis of [F] 

begin 1 

choose q E F with the lowest ld{g) w.r.t. -< 2 

r:={(<?,W(ff),0)}; Q:=0; G := {g} 3 

for each f £ F\{g} do 4 

Q:=Qu{(/,M(/),0)} 5 

repeat 6 

h:=0 7 

while Q 7^ and h = do 8 

choose g in {g, 9,P) £ Q with the lowest ld{g) w.r.t. -< 9 

Q:=Q\{{g,e,P)} 10 

if Cnterzon(.g, 6*, T) is false then h := NFL{g,G) 11 

if /i^O then G:=GU{h} 12 

if ld{h) = ld{g) then T ■.= T\j{{h,e,Pf\NML{h,G))} 13 

else T := TU {(/i,W(/i),0)} 14 

for each f in {f,-d, S) e T s.t. ld{f) y ld{h) do 15 

T:=T\{(/,^9,5)}; Q := Q U {(/, ^9, 5)}; G:=G\{/} 16 

for each {f,^,S) eT do 17 

T := T \ {(/, ^, 5)} U {(/, t?, 5 n NMlU, G))} 18 
while exist {g, e,P) eT and a; e NML{g, G)\P and, if g ^ 0, 19 

s.t. Zd(5^ • 5) ^ W(/) for ah / in (/, ^,S) £Q do 20 

choose such {g, 9,P),x with the lowest ld{dx • g) w.r.t. -<c 21 

T:=T\{(g,0,P)}U{(g,0,PU{x})} 22 

if CriieWon(5a; • g, 61, T) is false then h := NFL{d^ ■ g,G) 23 

if /i 9^ then G := G U [h] 24 

a ld{h) = ld{d^ ■ g) then T:=TU{(/i,6',0)} 25 

else T:=ru{(/i,W(/i),0)} 26 

for each f m {f.d.S) (^T with ld{f) ld{h) do 27 
T:=T\ {(/,,?, 5)}; Q := Q U {(/, ^}); G := G \ {/} 28 

for each (/,??, 5) e T do 29 

T := T \ {(/, ^, 5)} U {(/, ^,Sn NMLif, G))} 30 

until g ^ 31 

end 32 



Criterion{g, 9, T) is true if there is (/, '&,S) e T such that ld{g) = dLld{f) 
and lcm{6,'&) -<c ld{g). 
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Example 21. Q The well-known Janet example with three independent and 
one dependent variables (n = 3, m 1): 

( diiy ~ X2d33y = 0, 
\ d22y = . 

The above completion algorithm applied for Janet, Pommaret and Lexico- 
graphically induced divisions gives the following involutive bases, which coin- 
cide for both pure lexicographical and graded lexicographical main rankings 
compatible with ^ and which sorted in the descending lexicographical order: 



Grobner 
basis 


Involutive Bases 


Janet & Pommaret 


Lex. Induced 


diiy - X2d33y 


diiy - X2d33y 


diuy - dssy 


d22y 


di22y 


9ll333y 


9233 y 


9l2332/ 


9ll332/ 


933332/ 


9l3333?/ 


diisy- X2d333y 




d22y 


diiy - X2d33y 




d233y 


92232/ 




933332/ 


9222/ 






923332/ 






92332/ 






933332/ 



The first column contains the reduced differential Grobner basis, and Janet 
and Pommaret bases are identical for this example. 

6 Initial Value Problem 

The results of this section generalize to arbitrary L— involutive linear systems 
those obtained in Riquicr- Janet theory ||6|,|8| , ^ , for Janet and Thomas divi- 
sions, as well as in the formal theory [[l2yi3|] for Pommaret division, on posing 
an initial value problem providing uniqueness and existence of solutions. 

Definition 22. [U| If e ld{F) is a leading derivative in F C M, then d^O 
is called a principal derivative. A derivative which is not principal is called 
parametric. The monomial set W — {U^'Lj^Wj | Wj C M} associated by (0) 
with the set of parametric derivatives is called a complementary set of F. 

Proposition 23. Given a ranking if set F is a linear L— involutive basis 
of differential ideal [F] , then the sets of principal and parametric derivatives 
{complementary set) related to F depend only on [F] and -< and do not depend 
on the choice of involutive division L. 



Proof. It follows immediately from the fact that any involutive basis is a 
Grobner basis (Theorem . □ 
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Lemma 24. {decomposition lemma) Given a noetherian division L and L— 
involutive set _F C M, every subset Wj in the complementary monomial set of 
F related to j-th differential indeterminate yj {I < j < m) can be decomposed 
as a disjoint union 

Wj = U„ey^wL^, L„ C L{v, Uj U {«}), (5) 

where Uj is the L— involutive monomial set (not necessarily nonempty) asso- 
ciated with ldj{F), and Vj £ Wl is a finite subset. 

Proof. Let U and W he a pair of monomial sets associated with the prin- 
cipal and parametric derivatives of a differential indeterminate in F. The 
complementary set W can be written as a disjoint union [3^ 

W ^WoUWiU---UWd (6) 

where d is the dimension of monomial ideal ([/), Wq is a finite set, and every 
Wi (1 < « < d) is a finite disjoint unionQ 

Wr =Wr,UWr^U---U Wr^ (7) 

with 

Wr, = { Wr,<' • • • a;";J at e N , 1 < i < r } (1 < s < fc). (8) 

For every v e Wq we shall take = {1} in (||). Thus, for d = the 
decomposition (||) = Wq = Uv£Wo{v} holds trivially. If d > we consider 
the finite set 

V = WoUtl^s=l{'^rJ, (9) 

where monomials Wr^ generate Wr^ in accordance with (||) . 

We claim that elements in set (^), and the decompositions (||), (0) they 
determine can be written such that the union in = U^gywL^ with Ly C 
L{v, UL) {v}) is disjoint in accordance with (^. To prove the claim we define 
the degree q of set U a.s q — maxjdeg{u) \ u 6 U}, and choose all the 
monomials Wr^ generating Wr^ in such that deg{wr^) — q. Obviously 
this can always be done by appropriate choice of Wq. Let now Vi be the set 
Vi = Uf^i ^^=1 {wr^}, and let U he a finite L— autoreduced completion of 
UUV. The existence of U is guaranteed by noetherity of L. Now consider the 



set V = U nW 3 Vi. Its L— involutivity and property (d) of L in Definition 10 
imply 

{yweW\ Wo) {3veV) [we vL{v, T) C vL{v, U U {v}) ]. 

Thus, we obtain the desired decomposition W = Wo U^^y vL{v,U). Dis- 
jointedness of this union follows from that in (Q) and Definition || of L— 
autoreduction. This proves the claim and the lemma. □ 

^ The union in (^) considered in js^] is not necessarily disjoint. However, unions in 
(^) and (^) apparently can be rewritten as disjoint by appropriate choice of Wo 
and components of Wr- 
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Definition 25. Those elements Vj^ (parametric derivatives) which, in ac- 
cordance with (H), generate the whole complementary set W, will be called 
L— generators of the set. The multiplicative variables Xi satisfying Xi G Lj^, 
will be called [L—) multipliers of the generator Vj^ and the remaining vari- 
ables will be called its {L—)nonmultipliers. The whole set of i— generators of 
W will be denoted hy Gl, and in accordance with (||) 

GL^Uf^.V,. (10) 



For a non-noetherian division L a complementary set may not have a finite 
set of L— generators as the following example shows. 

Example 26. Let involutive division L be defined on M as follows. Vari- 
ables xi,...,Xn-i are separated into multiplicative and nonmultiplicative 
by Definition ^. Let the variable Xn be also separated by Definition ^ if 
degn{u) = and u ^ 1, whereas if degn{u) > or if u = 1, a;„ be nonmul- 
tiplicative for u. Then, the monomial set U = {x^, xiX2, X2} is L— involutive 
in K[xi, X2, X3]. Its complementary set has the infinite set of i— generators: 
GL = {l}U{x,}UfZ,{xl}. 



Remark 27. Decomposition (^) and the underlying L— generator set (£^) are 
not uniquely defined, and usually a more compact set G of L— generators (with 
less number of elements) than that constructed in the proof of Lemma G4| can 
be chosen. For example, for a Janet basis, Gp can always be chosen [pj as 
union (|l^) of sets Vj such that 

(VVj) {VveV,) [L, = Jiv,UjU{v}] (l<j<m), (11) 

where J stands for the Janet set of multiplicative monomials. Since for Uj, as 
it constructed in the proof, the inclusion Uj U {v} C Uj holds, the property 
(d) in Definition ^ imphes J{v,Uj) C J{v,Uj U {v}). Therefore, the set of 
Janet generators defined by (|l^) is a subset of that constructed in the proof 
of Lemma ^ 

For a Pommaret basis in the formal theory decomposition (||) is taken 
in the form 

W = Wa Ul^(,w I deq{v)=q} vP{v), (12) 

where P{v) denotes the set of Pommaret multiplicative monomials for v, 
and q, as in the proof, is the degree of the basis. The number of Pommaret 
generators in ( |l^ ) with i multipliers is called the ith Cartan characte^ (1 ^ 
t < n) of the basis and will be denoted by a^. 

^ Cartan introduced these numbers in his analysis of exterior PDEs Q and called 
them characters. 
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Example 28. The complementary set of the monomial ideal {U) for U = 
{xlx'i,xiX2,xixl} in Example |ll| is W = \J{x^'^ \ i G N}U{a;^24 \j,ke N}. 
Its most compact sets Qj and Qdl^^ together with their multipliers are: 



Janet division 


Lex. induced division 


Generator 


Multipliers 


Generator 


Multiphers 


1 


22,2:3 


1 


2:2, 23 


Xl 




21 


2:1 


x\ 


21 







We note that if the involutive bases Up^Udl^^ given in Example |ll| are 
sequentially enlarged with every single generator, then the sets of Janet mul- 
tipliers, in accordance with Remark ||, coincide with the sets of multiplicative 
variables 

Mj(l,;7pU{l})-{x2,23}, Mj(xi,C7pU{l}) = 0, Mj{xl,UpU{xl}) = {x,} 

whereas for lexicographically induced division, every set of multipliers is the 
proper subset of multiplicative variables 

Mdz^^A'^iUdz,,^ U {1}) = Mo^^A^uUDr^ex U {21}) = {21,22,2:3}. 

Theorem 29. {uniqueness theorem) Let J- be an L^involutive system of lin- 
ear PDEs for an orderly ranking. Then T has at most one solution satisfying 
the following initial conditions: the derivatives associated with L— generators 
of the complementary monomials are arbitrary functions of their multipli- 
ers at the fixed values of their nonmultipliers from coordinates of the initial 
point Xi = x° (1 < i < n)), whereas the generators without multipliers are 
considered to be arbitrary constants. 

Proof. Involutivity of J- with respect to an orderly ranking implies that the 
associated complementary monomial set contains all the monomials associ- 
ated with the parametric derivatives. This statement is an immediate con- 
sequence of the well-known fact ||3^ that for a graded monomial ordering 
the Hilbert function of a polynomial ideal is defined by the monomial ideal 
generated by the leading monomials of a Grobner basis of the polynomial 
ideal. 

Furthermore, by association (|l|), the decomposition (^) yields that every 
parametric derivative associated with a monomial in \ Wq is produced by 
differentiation of the uniquely defined parametric derivative (L— generator) 
with respect to its multipliers. Assigning the fixed values to all these paramet- 
ric derivatives is obviously equivalent to fixing some function of the multipli- 
ers. Therefore, given initial point Xi = x°, in addition to the set of arbitrary 
constants which associated with elements in Wq , all the parametric arbitrari- 
ness is determined by functions corresponding to the L— generators and which 
are arbitrary functions of the multipliers at the fixed values of nonmultipliers 
from coordinates of the initial point. □ 
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Remark 30. Different involutive divisions give obviously equivalent forms of 
initial value problem providing the uniqueness of solutions. However, given a 
system of PDEs with an infinite set of parametric derivatives, the writingf of 
such initial conditions in accordance with Theorem may be more compact 
for one division than for another. Wc demonstrate this fact by examples given 
below. 



Theorem 31. {existence theorem) Let J- he an L— involutive linear system 
for an orderly ranking, and let its coefficients he analytic functions in an 
initial point {xi = x°). Then T has precisely one solution which is analytic in 
this point if all arhitrary functions in the initial data specified in Theorem 
are analytic in their arguments taking values from coordinates of the initial 
point. 

Proof. This is identical to the existence proof in Riquier- Janet theory 
(see also 31 ). □ 



Example 32. The complementary monomial set for Janet system in Exam- 
ple ^ is finite and consists of 12 elements 



By Theorem 31 , its general solution depends on 12 arbitrary constants. 



Example 33. The system of the first order PDEs with four independent 
and one dependent variables {n = 4,,m— 1) and its completion to involution 
for Janet or Pommaret division for any ranking compatible with (|^) are given 

by 

d^y + x,d,y + y^O, ^,P-co„.pic«on (d,y + x,<hy + y = 0, 
d.y + x,d,y = 0, =^ ld,y + x,d,y^O, 

The parametric derivatives d\y {i G N) have the only Janet generator y <^=> 1 
with the only multiplier x^. Hence, the initial data providing the unique ana- 
lytic solution are y\xi=x°,x2=x°,x3=x° = (f'ixi) with arbitrary function (j){x4^), 
analytic at x^ — x^. The system can explicitly be integrated and its general 
solution is 

y = e"^^''~''^)(/)(a;3 +X4- X1X2 + X°X2 - 2:3). 

The Pommaret generators are y and d^y without multipliers and with mul- 
tiplier X4, respectively. This leads to the initial value problem 

y\xi—x°,X2—x'^,X3—x'^ — C, d4y\xi—x°,X2—x^,x:i—x'^ — ^^(^4) 

with arbitrary constant c and arbitrary function ^p. This shows that the 
Janet initial conditions are written in a more compact form than those of 
Pommaret. 
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Example 34- [ p7[ The well-known Lewy example with n = 3, m = 2 and 
r]i,r]2 e K 



This system is involutive for any of Janet, Pommaret or lexicographically 
induced divisions and the orderly ranking with diyj >~ d2yj >~ d^yj, yi >- y2- 
Janet generators are yi,y2- Each of them has multipliers X2,X3. This implies 
the initial data providing the uniqueness: yj\xi=x" — 4'jix2,X3) {j = 1,2) 
with arbitrary functions (j)j{x2, x^). Pommaret and lexicographically induced 
divisions lead to a less compact writing of these conditions. 

Remark 35. As shown by Lewy for Example ^ there exist the C°° 
functions f7i,»72 such that the system has no C°° (and even C^) solutions. 
Therefore, analyticity in the Theorem ^ statement can not be replaced by 
smoothness. 

We conclude this section with explicit formulae for the Hilbert function HF[p] 
and Hilbert polynomial HP[f] of differential ideal [F] represented by its linear 
involutive basis F . These formulae are valid for any involutive division and an 
orderly ranking. For ordinary differential ideals, that is, for the case of single 
differential indeterminate (to = 1), by association (^, they are the same as 
in commutative algebra [ p0|j22| . For partial differential case they involve the 
number m of differential indctcrminates 



i7P[,,(,s)^TO("+-^) E (14) 

J — 1 U j 



Here n is the number of independent variables, Uj is the monomial set asso- 
ciated with the set of leading derivatives ldj{F), and /i(w) is the number of 
multiplicative elements of u. 

The first term in the right hand side of (^3|) is the total number of deriva- 
tives of order < s. The triple sum counts the number of principal derivatives 
among them in accordance with Definition |l^ which says that any princi- 
pal derivative is uniquely obtained by the multiplicative prolongation of one 
of the leading derivatives in F . Thus, ( [l3| ) gives the number of parametric 
derivatives of order < s, and for s large enough it becomes polynomial (|lj). 

In the formal theory [ p"2|JT^ the Janet formula is used: 



{ 



-2x3 922/1 -dzy2 - 2xi 522/2 = ?/i(a;i, a;2, xa), 
diy2 + 2xi diyi + d^yi - 2x3^22/2 = f/2(a;i, a;2, 2:3). 
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Here the Hilbert polynomial ||] is written in terms of Cartan characters 
(see Remark^). Apparently, this is (14), rewritten for Pommaret division in 
terms of Cartan characters. 



7 Lie Symmetry Analysis of PDEs 

Lie symmetry methods and their computerization yield a powerful practical 
tool for analysis of nonlinear differential equations (see the review article [ p8| 
and references therein for more details). We present here the basic compu- 
tational formulae and demonstrate, by two simple examples with a single 
nonlinear evolution equation, application of the above described involutive 
methods to finding the classical infinitesimal symmetries. 
Given a finite system of polynomial-nonlinear PDEs 

fk{xi,yj,. . . ,daVj) = 0, {1 < i < n, I < j < m, 1 < k < r) (15) 

one looks for one-parameter infinitesimal transformations 

x.(A) = X. +6(x.,y,)A + 0(n (1 < , < 1 < . < (16) 



The conditions of invariance of (15) under transformations (|16|) are 

Z^"\fk{x^,yJ,...,^o.y,)\f^=o = 0, (l<fc,s<r) (17) 

= CA, + Vjdy, + Qndy,,, +■■■ + Cr.cdy,,^ , (18) 

where djy j d enoted by yj;i, etc.^ Functions Cj;... involved in the differential 
operator ( [l8[ ) are uniquely computed in terms of functions ^i,rij and their 
derivatives by means of the recurrence relations 

=D^{r|j) ~yr,qD^{£.q), 
Cj;ii...ip = ^ip(Ci;ii---4p-i) ~ yj;ii . . .ip_iqDi^{£_q) , 

where Di is the total derivative operator with respect to Xi 

Di = di + yj-idy^ + yj-ikdy..^ H 

The invariance conditions (^7|) produce the overdetermined system of lin- 
ear homogeneous PDEs in 5i,77j which is called the determining system. Its 
particular solution yields an infinitesimal operator of the symmetry group 

Z ^C^d^,+T]jdy^, (19) 

and the general solution yields all the infinitesimal operators. 

Given initial system ([l5|), integration of the determining system is gen- 
erally a bottleneck of the whole procedure of constructing these symmetry 
operators, and completion the system to involution is the most universal 
algorithmic method of its integration . 

^ In this section the summation over repeated indices is always assumed. 
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Example 36. p3[ Diffusion type equation yt + yy-x — tyxx = 0(n = 2,m=l). 
The symmetry operator iMm of the form 



Z^^idt + i2dx + vdv (20) 
satisfies the determining system 

dyyil = 0, dyyi2 = 0, t dyyT] ^ 2 t dxyS,2 " 2 ^ dy^2 = 0, 

dy^i = 0, 2 t^a^^,^ _ t^dxx^2 ~ytdx^2 + t dt^2 + 2/ a - * = 0, 

t dxxV - y dxV -dtr]^0, dxxii - yt dx^i +2t dx^2 - t dt^i - = 0, 

tdxyS.1 + dy£_2 = 0, dx^i 0. 

By choosing the orderly degree-reverse-lexicographical ranking with dy >- 
dx ^ dt, £,1 )^ ^2 ^ V and applying the completion algorithm of Sect. 5, we 
obtain the (Pommaret, Janet, lexicographically induced) involutive system 

dy^i^O, 9^6 = 0, 9,77 = 0, 9,6=0, 9,6-^6=0, 

dxV^O, 9*6-^^1=0, dt£2-V^0, dtV^O. 

The generators of parametric derivatives £,i,£,2,'ri have no multipliers. Hence, 
the general solution depends on three arbitrary constants ci, C2, C3, and it can 
easily be obtained by explicit integration of the involutive system 

6 = cii, £2 ^ cix + C2t + C3, ri = C2. 

Respectively, the Lie symmetry group is three-dimensional. Its symmetry 
operators Zi = tdt + xdx, Z2 = tdx + dy, — dx form the Lie algebra 
[Zi, Z2] = 0, [Z2, Za] = 0, [Z1J3] = -Z3. 

Example 37. The Harry Dym equation dty — y^dxxxy = {n = 2, m = 1) 
which was already used in as an illustrative example. The symmetry 
operator in the form ( ^ ) is now determined by the system 

dy^l = 0, dx^l = 0, dy£,2 = 0, dyyT] = 0, 

dxyV - dxx£,2 = 0, dtrj - y^dxxxV = 0, 

3 y^dxxyV + dt(,2 - y^dxxx£.2 =0, y dt£.i - 3 y dx£.2 +3t] = 0. 

Its Janet and Pommaret involutive form for the same ranking as in the pre- 
vious example is 

OxxV = 0, dxtV = 0, dyf] - ^1] = 0, dtrj = 0, = 0, 

9.6 - ^ 9t6 - - ^7 = 0, 9t6 = 0, 9„ei = 0, dy^i = 0, 9,6 = 0. 
3 y 
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There are five generators of parametric derivatives ^i,dt^i,S,2,'ri,dx'fj wliich 
have no multiphers that implies the five-dimensional Lie symmetry group. 
The involutive determining system in this example is also easy to integrate: 

^i=ci+C2t, 6 = C3 + C4X + csx^, ri = {c4 - -C2 + 2c5x)y. 
This gives the Lie symmetry operators 

Zi=dt, Z2^tdt-^ydy, Zs^d^, Z^^xd^+ydy, Z5 ^ x'^dx + 2xy dy 
with the following nonzero commutators of the symmetry algebra 

[Zi,Z2]=Zi, [Z3,Z4^Zi, [Z3,Z5]=2Z4, [Zi,Z5]^Z5. 



8 Conclusion 

Most of the above presented definitions, statements and constructive meth- 
ods can be extended to finite sets of differential polynomials in M which, 
given a ranking, are linear with respect to their highest rank (principal) 
derivatives. In Riquier- Janet theory the corresponding systems of PDEs are 
called orthonomic. Their completion to involution, for any constructive and 
noetherian division, could be done much like linear systems. The essential ob- 
struction here is a non-orthonomic integrability condition. Moreover, even if 
such an integrability condition is explicitly solvable with respect to its princi- 
pal derivative, then this leads to non-polynomial orthonomicity, and, thereby, 
to difficulty in the use of constructive methods of differential and commuta- 
tive algebra. In the latter case some geometric features of the formal theory 
may be useful for computational purposes . 

However, given an orthonomic system of polynomial PDEs and an involu- 
tive division L, one can always verify if it is i— involutive. Analytic involutive 
orthonomic systems admit posing an initial value problem providing the ex- 
istence and uniqueness of solution. One can, hence, determine arbitrariness 
in the general solution as it is done in Sect. 6 for linear systems. In particu- 
lar, the compact general formulae ( p^ ) and (|l^) for the Hilbert function and 
Hilbert polynomial are also valid for involutive orthonomic equations. 

We are going to implement the completion algorithm MinimalLinear- 
InvolutiveBasis ( Sect. 5 ) after examination and optimization of its poly- 
nomial analogue pl| , |23[ | . Though its implementation in Reduce for Pommarct 
division has already shown its efficiency, the differential case needs more 
careful analysis of implementation and optimization issues to be applicable 
to PDEs of practical interest. Thus, in Lie symmetry analysis of relatively 
small systems it is easy to obtain determining systems of many hundreds 



21 



equations. Currently, the most efficient completion algorithm for linear sys- 
tems implemented in some packages for Lie symmetry analysis [ p8| is that 
of paper . Its underlying implementations allow to treate hundreds deter- 
mining equations (cf. [p5[). As for significantly larger determining systems, 
they are hardly tractable by the present day computer algebra tools, whereas 
there are practical needs in it. In gas dynamics, for instance, the group clas- 
sification of the system of five second order PDEs describing a viscous heat 
conducting gas and involving five dependent and four independent variables 
(three spatial and one temporal) [^|, leads to the determining system con- 
taining more than 200 000 equations. 

In our intention to extract, in the process of implementation, the maximal 
possible efficiency from the algorithms proposed, we hope, first of all, to detect 
(heuristically) the most optimal choice of involutive division. As the first step 
in this direction an implementation of the monomial completion for different 
divisions has been done in Mathematica and for Janet division in C IP4|. 
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